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In our previously published papers, it was argued that a massive non-Abelian 
gauge field theory in which all gauge fields have the same mass can well be set 
up on the gauge-invariance principle. The quantization of the fields was per- 
formed by different methods. In this paper, It is proved that the quantum theory 
is invariant with respect to a kind of BRST-transformations. From the BRST- 
invariance of the theory , the Ward-Takahashi identities satisfied by the generat- 
ing functionals of full Green's functions, connected Green's functions and proper 
vertex functions are successively derived. As an application of the above Ward- 
Takahashi identity, the Ward-Takahashi identity obeyed by the massive gauge 
boson propagator is derived and the renormalization of the propagator is dis- 
cussed. Furthermore, based on the Ward-Takahashi identity, it is exactly proved 
' that the S-matrix elements given by the quantum theory are gauge-independent 

, and hence unitary. 

■ PACS: 11.15-q, 12.38-t 

m ! 

> ■ 1. INTRODUCTION 

(N ; 

^| ' In our preprevious papers' 1-4 ', it has been shown that some massive gauge field theory in 

■ which the masses of all gauge fields are the same may really be set up on the gauge-invariance 
principle without the help of the Higgs mechanism. The essential points of achieving this con- 
clusion are the following. (1) The massive gauge field must be viewed as a constrained system in 
the whole space of vector potential. Therefore, the Lorentz condition, as a necessary constraint, 
must be introduced from the onset and imposed on the massive Yang- Mills Lagrangian so as to 

+f* i restrict the unphysical degrees of freedom involved in the Lagrangian; (2) The gauge-invariance 
Q_ij of gauge field dynamics should be more generally required to the action of the field other than 
the Lagrangian because the action is of more fundamental dynamical meaning. In particular, 
' the gauge-invariance for the constrained system should be required to the action written in the 
;> physical subspace defined by the Lorentz condition in which the fields exist and move only; (3) 
In the physical subspace , only the infinitesimal gauge transformations are possible to exist and 
necessary to be considered in inspection of whether the theory is gauge-invariant or not; (4) 
To construct a correct gauge field theory, the residual gauge degrees of freedom existing in the 
physical subspace must be eliminated by the constraint condition on the gauge group. This 
constraint condition may be determined by requiring the action to be gauge-invariant. Thus, 
the theory is set up from beginning to end on the gauge-invariance principle. These points 
are important to build up a correct quantum massive non-Abelian gauge field theory. Such a 
theory, as will be proved, is renormalizable and unitary. 

In Refs. [1-3], the quantum theory of the massive non-Abelian gauge fields without Higgs 
mechanism was established by different methods of quantization. In this paper, it will be 
shown that the quantum theory has an important property that the effective action appearing 
in the generating functional of Green's functions is invariant with respect to a kind of BRST- 
transformations' 5 '. From the BRST-symmetry, we will derive various Ward-Takahashi (W-T) 
identities' 6-12 ! satisfied by the generating functionals of Green's functions and proper vertices. 
These W-T identities are of special importance in proofs of unitarity and renormalizability of 
the theory. In this paper, we confine ourselves to prove that the S-matrix elements evaluated 
from the theory is independent of the gauge parameter, that is to say, the gauge-dependent 
unphysical poles appearing in the gauge boson propagator and the ghost particle propagator 
do not contribute to the S-matrix elements. Therefore, the unitarity of the theory is ensured. 
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The arrangement of this paper is as follows. In section 2, we will derive the BRST- 
transformations under which the effective action of the massive non-Abclian gauge field theory 
is invariant. In doing this, we extend our discussion by including fermions. In section 3, we will 
derive the W-T identities satisfied by various generating functionals. In section 4, to illustrate 
applications of the above W-T identity, the W-T identity obeyed by the massive gluon prop- 
agator will be derived and the rcnormalization of the propagator will be discussed. In section 
5, by virtue of the W-T identity, we will prove that the S-matrix elements calculated from the 
massive gauge field theory are gauge-independent and hence unitary. The last section is used 
to make some remarks on the nilpotency problem of the BRST-transformations and the BRST- 
invariace of the external source terms introduced to the generating functional. In Appendix, 
the W-T identities used to prove the unitarity will be given by an alternative derivation. 



2. BRST- TRANSFORMATION 



In the previous paper , we mainly discussed the gauge fields themselves without concerning 
fermion fields. For the gauge fields, in order to guarantee the mass term in the action to be 
gauge- invariant, the masses of all gauge fields are taken to be the same. If fermions are included, 
as pointed out in Rcfs.[l-3], the QCD with massive gluons fulfils this requirement because all 
the gluons can be considered to have the same mass m. The SU(3)-symmetric action of the 
QCD with massive gluons is given by the following Lagrangian^ 1-4 ! 

C = ^(d, igT a AV - M}yj - -F^F^ + \m 2 A a » A^ (2.1) 

where i[)(x) denotes the quark field function , 4>{x) is its Dirac-conjugate, T a = X a /2 are the 
color matrices and M is the quark mass. The above Lagrangian is constrained by the Lorentz 
condition 



d^A* = (2.2) 

Under this condition, as was proved in Refs. [1-3], the action given by the Lagrangian in Eq. 
(2.1) is invariant with respect to the following gauge transformations: 

5A- = iDfC h 

5ip_{x) = ig£T a C a (x)i{>(x) (2.3) 
5ip{x) = ig&(x)T a C a (x) 

where we have set the parametric functions of the gauge group a (x) — ^C a {x) in which £ is 
an infinitesimal Grassmann number and C a {x) are the ghost field functions. According to the 
result given in paper I, the quantum theory built up from the Lagrangian in Eq. (2.1) and the 
constraint condition in Eq. (2.2) is described by the following generating functional of Green's 
functions! 5-9 ] 



Z[J^K a ,K a ,fj,T}] ^^jV(Al 1 C a 1 C a ^,,p)exp{tS 
+i J d 4 x[J a ^Al + K a C a + C a K a + rpr) + #]} 

where J aAt , K a , K a , 77 and fj are the external sources and 



(2.4) 



S = J d 4 x{$[i^(d^ - igT a Al) - M]V - \F a ^F« u + \vr?A a »Al 

-h^ A V) 2 + c a d»{vfc b )} ' 1 ' 

is the effective action in which 

Vf{x)=>^d* + Df{x) (2.6) 

here /i 2 = am 2 and 

Df = 5 ab d„ - gf abc Al (2.7) 
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is the covariant derivative. Similar to the massless gauge theory, for the massive gauge the- 
ory, there are a set of BRST-transformations including the infinitesimal gauge transformations 
shown in Eq. (2.3) and the transformations for the ghost fields under which the effective action 
is invariant. The transformations for the ghost fields may be found from the stationary condi- 
tion of the effective action under the BRST-transformations. By applying the transformations 
in Eq. (2.3) to the action in Eq. (2.5), one can derive 

SS = J d 4 x{[5C a - ^d v A a v ]d^{V ab C b ) + C a d»5{VfC h )} = (2.8) 

This expression suggests that if we set 

SC a = -d v A a v (2.9) 
a 

and 

& t 5(pfC b )=Q (2.10) 

The action will be invariant. Eq. (2.9) gives the transformation law of the ghost field variable 
C a (x) which is the same as the one in the massless gauge field theory. From Eq. (2.10), we may 
derive a transformation law of the ghost variables C a (x). Noticing the relation in Eq. (2.6), we 
can write 

S(Vf(x)C b {x)) = ^d*5C a (x) + 5{Df{x)C\x)) (2.11) 

In the massless gauge theory, it has been proved that' 8-11 ! 

S(Df{x)C\x)) = Df(x)[5C\x) + ^gf bcd C c (x)C d (x)} (2.12) 

With this result, Eq. (2.11) can be written as 

5{Vf{x)C\x)) = Vf{x)5C\x) - Df{x)8C b {x) (2.13) 

where 

SCS(x) = -^f abc C b (x)C c (x) (2.14) 

On substituting Eq. (2.13) into Eq. (2.10), we have 

M ab (x)5C b (x) = Mg b (x)5C%(x) (2.15) 

where we have defined 

M ab (x) = d%V ab {x) = S ab (D x + M 2 ) - ff/ o6c A= (s)S^ (2.16) 

and 

M afc (x) = d»D ab {x) = M ab {x) - fi 2 6 ab (2.17) 

It is noted that the operator in Eq. (2.16 ) is just the operator appearing in the equation of 
motion for the ghost field C a (x) 

Qi*(priC b )=0 (2.18) 

(see Eq. (4.10) in paper I). Corresponding to this equation of motion, we may write an equation 
satisfied by the Green's function A ab (x — y) 

M ac (x)A cb (x - y) = S ab S 4 (x - y) (2.19) 
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The function A ab (x — y) is nothing but the exact propagator of the ghost field which is inverse 
of the operator M ab (x). In the light of Eq. (2.19) and noticing Eq. (2.17), we may solve out 
the 5C a {x) from Eq. (2.15) 

5C a (x) = {M-^MoSCoYix) = {M~\M - ^)5C } a {x) 

= SCS(x) - m 2 / d*yA ab (x - y)SC b (y) l ^ Uj 

This just is the transformation law for the ghost variables C a (x). When the mass tends to zero, 
Eq. (2.20) immediately goes over to the corresponding transformation given in the massless 
gauge field theory. It is interesting that in the Landau gauge (a = 0), due to \i = 0, the above 
transformation also reduces to the form as given in the massless theory. This result is natural 
since in the Landau gauge, the gauge field mass term in the action is gauge-invariant. However, 
in general gauges, the mass term is no longer gauge- invariant. In this case, to maintain the 
action to be gauge- invariant, it is necessary to give the ghost field a mass /x so as to counteract 
the gauge-non-invariance of the gauge field mass term. As a result, in the transformation given 
in Eq. (2.20) appears a term proportional to [i 2 . 



3. WARD-TAKAHASHI IDENTITIES 



This section serves to derive the W-T identities for the quantum massive non-Abclian gauge 
field theory established in paper I and represented in Eqs. (2.4)-(2.7) on the basis of the BRST- 
symmetry of the theory. Since derivations of the W-T identities for the QCD with massive 
gluons are much similar to those for the QCD with massless gluons, we only need here to give 
a brief description of the derivations. When we make the BRST-transformations shown in Eqs. 
(2.3), (2.9) and (2.20) to the generating functional in Eq.(2.4) and consider the invariance of 
the generating functional, the action and the integration measure under the transformations 
(the invariance of the integration measure is easy to check), we obtain an identity such that 

i / V(A«, C a , C a , V) / d 4 x{ J a »(x)5A*{x) + 5C a {x)K a {x) + K a (x)6C a (x) 

+fj(x)6ip(x) + S^(x)r)(x)}e tS+EST (3.1) 
= 

where EST is an abbreviation of the external source terms appearing in Eq. (2.4). The 
Grassmann number £ contained in the BRST-transformations in Eq. (3.1) may be eliminated 
by performing a partial differentiation of Eq. (3.1) with respect to £. As a result, we get a W-T 
identity as follows 

^ / V(A°,C a ,C a ,ij>,1>) J d 4 x{J a »{x)AA<^{x) + AC a (x)K a (x) - K a (x)AC a (x) 

-f)(x)AtP(x) + A^(x)r](x)}e iS + EST (3.2) 
= 



where 

AAl(x) = Df{x)C b {x) 
AC a (x) = ±<9^(x) 
AC a (x) = J d 4 y[<5 Qb <5 4 (x - y) - ^A ab (x - y)]AC b {y) . , 

AC b (y) = -\gf hcd C c (y)C d {y) ^ 
Af{x) = igT a C a (x)^(x) 
Aip(x) = igij(x)T a C a {x) 

These functions defined above are finite. Each of them differs from the corresponding BRST- 
transformation written in Eqs. (2.3), (2.9) and (2.20) by an infinitesimal Grassmann parameter 

t 

In order to represent the composite field functions AA 1 ^ , AC a , Aip and Aip in Eq. (3.2) 
in terms of differentials of the functional Z with respect to external sources, we may, as 
usual, construct a generalized generating functional by introducing new external sources (called 
BRST-sources later on) into the generating functional written in Eq. (2.4), as shown in the 
following^ 11 ! 
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z[j-,K a ,K a ,n,mu a ^v a xx] 

= ±J V_[Af l ,C a ,C a ,^,^}exp{iS + i J <i 4 x[u<^AA£ +y a AC a (3.4) 
+AiP( + CAV> + J ati Al + K a C a + C a K a + r)i}) + V"?]} 

where w aM , v a , c and <; arc the sources which belong to the corresponding functions AA 1 ^ AC a , 
A<J> and A\& respectively. Obviously, the w aM and AA® are anticommuting quantities, while, 
the v a , (, C, AC a , Af/S and Atj) are commuting ones. We may start from the above generating 
functional to re-derive the W-T identity. In order that the identity thus derived is identical 
to that as given in Eq. (3.2), it is necessary to require the BRST-source terms UiA&i where 
m = u a ^, v a , C or C and A$ 4 = AA a ^ AC a , A* or AW to be invariant under the BRST- 
transformations. How to ensure the BRST-invariance of the source terms? For illustration, let 
us introduce the source terms in such a fashion 

/ d 4 x[u a ^SA^ + v a 5C a + (Sip + frip(} , 3 5) 

= J d A x[u a ^AAl + v a AC a + (Atp + A^C] 

where 

u " = « ^, v a = v a £, ? = ? = (3.6) 

These external sources are defined by including the Grassmann number £ and hence products 
of them with £ vanish. This suggests that we may generally define the sources by the following 
condition 



Ui£ = (3.7) 

Considering that under the BRST-transformation, the variation of the composite field functions 
given in the general gauges can be represented in the form <5A$i = £<I>i where are functions 
without including the parameter £, clearly, the definition in Eq. (3.7) for the sources would 
guarantee the BRST- invariancc of the BRST-source terms. When the BRST-transformations 
in Eqs. (2.3), (2.9) and (2.20) are made to the generating functional in Eq. (3.4), due to the 
definition in Eq. (3.7) for the sources, we have UiSA^i — which means that the BRST-source 
terms give a vanishing contribution to the identity in Eq. (3.1). Therefore, we still obtain 
the identity as shown in Eq. (3.1) except that the external source terms is now extended to 
include the BRST-cxtcrnal source terms. This fact indicates that we may directly insert the 
BRST-source terms into the exponent in Eq. (3.1) without changing the identity itself. When 
performing a partial differentiation of the identity with respect to £, we obtain a W-T identity 
which is the same as written in Eq. (3.2) except that the BRST-source terms are now included 
in the identity. Therefore, Eq. (3.2) may be expressed as 

+»K*)*fe + h Ka (*WiAaWZ> ■ ■ ■ > C] (3-8) 

= 



This is the W-T identity satisfied by the generating functional of full Green functions. 
On substituting in Eq. (3.8) the relation 6-8 

Z = e tW (3.9) 

where W denotes the generating functional of connected Green's functions, one may obtain a 
W-T identity expressed by the functional W 

+ L K a [x)di _A_ ]W[J a i ... X] (3 . 10) 

= 



From this identity, one may get another W-T identity satisfied by the generating functional Y 
of proper (one-particle-irreducible) vertex functions. The functional T is usually defined by the 
following Legendre transformation! 8-11 ! 
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T[A a », C a , C\ i>, i>- u^v\ C, C] = W[J*, K a ,K a , fj, V ; u», v a , (, Q 
- J d A x[J*A a ^ + K a C a + C a K a + fjip + V^] 

where A", C a , C a , -tjj and ^ arc the field variables defined by the following functional 
derivatives' 5-8 ] 



Aa( \ _ sw r a <T\ — sw r a (T~\ — sw 

From Eq.(3.11), it is not difficult to get the inverse transformations' 5-8 ' 

J atJ -(x) = —g^ x y K a (x) — s J^ x yK a (x) = - s JF (x y 

It is obvious that 

SW__ST_ W_ST_ SW _ ST SW _ 5T 

Employing Eqs. (3.13) and (3.14), the W-T identity in Eq. (3.10) will be written as 

sr sr i sr sr , <5r sr 



(3.12) 



(3.13) 



(3.14) 



fd 4 x{ 



SA° (x) J^JF) 8C a (x) 8v{x) S^{x) 8((x) 

_sr 

8C°(x) 



+ (3-15) 







This is the W-T identity satisfied by the generating functional of proper vertex functions. 

The above identity may be represented in another form with the aid of the so-called ghost 
equation of motion. The ghost equation may easily be derived by firstly making the translation 
transformation: C a — ► C a + X a in Eq.(2.4) where A a is an arbitrary Grassmann variable, then 
differentiating Eq. (2.4) with respect to the A a and finally setting A = 0. The result is I 8-11 ! 

1 J D{A% C\ C a >, ^){K a (x) + d%(Vf(x)C b (x))}e iS+EST = (3.16) 

When we use the generating functional defined in Eq. (3.4) and notice the relation in Eq. (2.6), 
the above equation may be represented as' 8-11 ' 

[K a (x) - id£ -. \ , - ^/i 2 J Z[J° • • • , C] = (3.17) 

On substituting the relation in Eq. (3.9) into the above equation, we may write a ghost equation 
satisfied by the functional W such that 

K a {x)+d^ 5W - +M % r > VK ; =0 (3-18) 
v ' x 5u a ^{x) ^ 8K a (x) y ' 

From this equation, the ghost equation obeyed by the functional T is easy to be derived by 
virtue of Eqs. (3.12) - (3.14)! 8 - 11 ] 

d x X \ -H 2 C a (x)=Q (3.19) 



5C a {x) x Su a ^(x) 

Upon applying the above equation to the last term in Eq. (3.15). the identity in Eq. (3.15) 
will be rewritten as 

r,j4„r <sr sr , sr sr , sr sr , sr sr 

j a x\ 5A a Sua „ T SC a Sv a 8ip SQ Sip <5C 

+m 2 d v A a v C a - IdvpAZ-^} (3-20) 
= 
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Now, let us define a new functional f in such a manner 



2a 

From this definition, it follows that 



r =r / ,v.,io»A%y (s.21) 



1)1 >u +-d»d u Al (3.22) 



5A» 5A» a 

When inserting Eq. (3.21) into Eq. (3.20) and considering the relation in Eq. (3.22), we arrive 
at 



fd*x{ 



4„r sr 6r i sr sr , sr sr 



g A a Su a» I SC a Sv a ' Sip SC 

.rfsr +m *ff'AZC a } ( 3 - 23 ) 
= 



The ghost equation represented through the functional f is of the same form as Eq. (3.19) 

' d x l x -H 2 C a (x) =0 (3.24) 



SC a (x) x 5u a »(x) 

In the Landau gauge, since /x = and 9"j4" = 0, Eqs. (3.23) and (3.24) respectively reduce to 

6-8 

f 4 <5f $f if <5f Sfsf 5TST 

and 

$--8*^=0 (3.26) 

These equations formally are the same as those for the massless gauge field theory. 

Now, we would like to give another form of the W-T identity. The ghost equation (3.16) 
suggests that the first external source term in Eq.(3.5) which appearing in the generating 
functional in Eq.(3.4) may be replaced by 

J d i xu a ^{x)5Al{x) = J d 4 xu a "{x)AA^{x) (3.27) 

where 

SA^x) = iVf{x)C\x) = ZAA1(x) (3.28) 



with 

□ 



AA^(x) = AA^ + ^-d;C a (x) (3.29) 



In this case, from the above relations, we see, the W-T identity in Eq. (3.8) will be rewritten 



as 



^ja, {x) ^ d ,^ + l K a {x)d ,_^ }Z[ ja r .. X]=Q 

and the ghost equation in Eq. (3.17) becomes 

{^>)-^^)W^^]=0 (3.31) 
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Repeating the derivations described In Eqs .(3.9)-(3.15) and (3.17)-(3.24), one may obtain from 
Eqs. (3.30) and (3.31) the identities expressed by the functional T 



fj4„r sr dr , <5r sr , sr sr 

J U "H SA^{x) Su a *(x) "r" c5C a (iE) Sv a (x) ~r j/^xj <5C(x) 

+4fe)7§y - 3%£^0r) Wd<M«(z)C a (z)} (3-32) 
= 



(3.33) 



<5C°(x) x 6u a »(x) 

In comparison of Eqs. (3.32) and (3.33) with Eqs. (3.23) and (3.24), we see, the advantage of 
using 5A^ to define the external source is that the ghost equation (3.33) becomes homogeneous. 
However, the price paying for this advantage is the increase of an inhomogeneous term (the fifth 
term) in Eq. (3.32). In the Landau gauge and in the zero-mass limit, Eqs. (3.32) and (3.33) 
still reduce to the homogeneous equations (3.25) and (3.26). 

From the W-T identities formulated above, we may derive various W-T identities obeyed 
by Green's functions and vertices, as will be illustrated later on. Particularly, these identities 
provide a firm basis for the proof of renormalizability and unitarity problems of the quantum 
massive gauge field theory as will be discussed in this paper and the next paper. 



4. PROPAGATORS 

In this section, as an application of the W-T identities derived in the preceding section, we 
have a particular interest in deriving the W-T identities satisfied by the massive gluon and 
ghost particle propagators and then discussing their renormalization. To derive the mentioned 
W-T identities, it is appropriate to start from the W-T identity in Eq. (3.30) and the ghost 
equation in Eq. (3.31). While, we would rather here to use the corresponding identities shown 
in Eqs. (3.8) and (3.17). Let us perform differentiations of the identities represented in Eqs. 
(3.8) and (3.17) with respect to the external sources K a {x) and K b {y) respectively and then 
set all the sources except for the source J?(x) to be zero. In this way, we obtain the following 
identities 



1 ™ 5Z \ J \ f ,4 rb», \ S 2 Z\J,K,u] 



and 



■ os 5 2 Z[J.K.u] I, . - 2 S 2 Z[J,K,K] , _ 

i( > Su^(x)SK»(y) \K=u=0 f «M SK a (x)SK b (y) \K=K=0 2 ) 

+S ab S 4 (x - y)Z[J] = 

Furthermore, on differentiating Eq. (4.1) with respect to J„(y) and then letting the source J 
vanish, we may get an identity which is, in operator representation, of the form I 8-11 ! 

h% < 0+|T[i«(z)it(y)]|0- >=< + \T*[C"(x)Dt d (y)C d (y)}\0- > (4.3) 

where A®(x), C a (x) and C a (x) stand for the gauge field and ghost field operators and T* 
symbolizes the covariant time-ordering product. When the source J is set to vanish, Eq. (4.2) 
will give such an equation^ 11 ! 

id" < 0+\T*{C"(x)Dt d (y)C d (y)}\0- > 
+i f i 2 < 0+\T[C a (x)C b (yW~ >= S ab S 4 (x -y) 

Upon inserting Eq. (4.3) into Eq. (4.4), we have 

&WD$(x -y)- a^A ab (x - y) = -aS ab S 4 (x - y) (4.5) 
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where 

l D^( X -y)=<0+\T{A;( X )A b u (y)}\0- > (4.6) 

which is the familiar full gluon propagator and 

iA ab (x - y) =< 0+\T{C a (x)6 b (y)}\0- > (4.7) 

which is the full ghost particle propagator. Eq. (4.5) just is the W-T identity respected by the 
gluon propagator which establishes a relation between the longitudinal part of gluon propagator 
and the ghost particle propagator. Particularly, in the Landau gauge, Eq. (4.5) reduces to the 
form which exhibits the transversity of the gluon propagator. By the Fourier transformation, 
Eq. (4.5) will be converted to the form given in the momentum space as follows 

k»k v D ab v {k) - an 2 A ab (k) = -a8 ab (4.8) 

The ghost particle propagator may be determined by the ghost equation shown in Eq. (4.4). 
However, we would rather here to derive its expression from the Dyson-Schwinger equation! 13 ! 
satisfied by the propagator which may be established by the perturbation method. 

A ab (k) = Af(k) + A* a '(k)n a ' b '(k)A b ' b (k) (4.9) 

where 

— if) ab 

l Af(k) = iS" b A (k) = k2 _^ 2 + i£ (4-10) 

is the free ghost particle propagator obtained in paper I and —iQ ab (k) = —i5 ab Q(k) denotes 
the proper self-energy operator of ghost particle. From Eq. (4.9), it is easy to solve that 

— it) ab 

iA ab {k) = ; (4.11) 

k 2 [l + Q(k 2 )] -^+ie 

where the self-energy has properly been expressed as 

n(k) = k 2 Cl(k 2 ) (4.12) 

Similarly, we may write a Dyson-Schwinger equation for the gluon propagator by the pertur- 
bation procedure' 13 ' 

D„Ak) = D%(k) + DP x (k)n x P(k)D pv (k) (4.13) 
where the color indices are suppressed for simplicity and 

»I><°> o6 (Jfe) = iS ab D^(k) = -i5 ab [ 9 ^ ~ \ kv/k2 + 1 f» k : /k2 . ] (4.14) 
^ w ^ w l k 2 -m 2 + ie k 2 -y? + ie 

is the free gluon propagator as derived in paper I and — iTl ab v {k) = — iS ab U^ u (k) stands for 
the gluon proper self-energy operator . Let us decompose the propagator and the self-energy 
operator into transverse and longitudinal parts: 

D^{k)=D^{k)+D^{k),W v {k)=T^ v {k)+Tl^{k) (4.15) 

where 

DPW = Or - !^f)D T (k 2 ), Df(fc) = ^D L (k 2 ), 



T 



Considering these decompositions and the orthogonality between the transverse and longitudi- 
nal parts, Eq. (4.13) will be split into two equations 
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D Tl „,{k) = D° T ^(k) + D^ x (k)U x /(k)D Tp ,(k) (4.17) 

and 

D Lllv {k) = D° L ^(k) + D% x (k)IL x L "(k)D Lpu (k) (4.18) 
Solving the equations (4.17) and (4.18), one can get 

iD ab (k) - -i5 ab l 9lLV - k ^ ku l k2 + ak^K/k 2 

tu^KK) i k 2 +U T (k 2 )-m 2 +ze + fc 2 + aII L (fc 2 )- M 2 W' { ' 

With setting 

H T (k 2 ) = fc 2 IIi(A; 2 ) + m 2 Ii 2 {k 2 ) (4.20) 
which follows from thr Lorentz-covariance of the operator lTr(fc 2 ) and 

aU L (k 2 ) = k 2 U L (k 2 ), (4.21) 

Eq. (4.19) will be rewitten as 



l ab (k) = -?7> afc i k^kyjk 2 ak^k v /k 2 

M ' t fc 2 [i + n 1 (fc 2 )]-m 2 [i-n 2 (P)]-+ze + fc 2[i + fi i (fc2 ) ]_ A ,2 + ze - 



We would like to note that the expressions given in Eqs. (4.12), (4.20) and (4.21) can be 
verified by practical calculations and are important for renormalizations of the propogators and 
the gluon mass. 

Substitution of Eqs. (4.11) and (4.22) into Eq. (4.8) yields 

n L ( fc 2 } = M 2 »( fc2 ) (4.23) 
k^[l + Q(k 2 )} 

From this relation, we see, either in the Landau gauge or in the zero-mass limit, the IIl(A: 2 ) 
vanishes. 

Now let us discuss the renormalization. The function Cl(k 2 ) in Eq. (4.11), the functions 
Hi(k 2 ), II 2 (fc 2 ) and I1l(A: 2 ) in Eq. (4.22) are generally divergent in higher order perturba- 
tive calculations. According to the conventional procedure of renormalization, the divergences 
included in the functions Cl(k 2 ), IIi(£; 2 ), n 2 (fc 2 ) and Il^fc 2 ) may be subtracted at a renormal- 
ization point, say, k 2 = v 2 . Thus, we can write! 5-9 ! 

n(k 2 ) = n(v 2 ) + mk 2 ), n 1 (fc 2 ) = n 1 (^) + nf(fc 2 ), 
n 2 (fc 2 ) - rr> 2 ) + n^(fc 2 ), n L (fc 2 ) = ii L ( v 2 ) + ii%(k 2 ) 1 ' > 

where (l{v 2 ), IIi(^ 2 ), n 2 (^ 2 ), IIl(^ 2 ) and n c (k 2 ), Il^fc 2 ), Il^fc 2 ), fl c L (k 2 ) are respectively the 
divergent parts and the finite parts of the functions Q(k 2 ), IIi(fc 2 ), II 2 (fc 2 ) and ili(fc 2 ). The 
divergent parts can be absorbed in the following renormalization constants defined byl 5 ~ 9 l 



z- 1 = i + n(v 2 ), z- 1 = 1 + i l> 2 ), z'- 1 = i + n L (t/ 2 ), 
z m ^ = y/z 3 [i - n 2 (^ 2 )] = vTi - ni(^ 2 )][i - n 2 (^ 2 )] 



(4.25) 



where Z 3 and Z 3 are the renormalization constants of gluon and ghost particle propagators 
respectively, Z 3 is the additional renormalization constant of the longitudinal part of gluon 
propagator and Z m is the renormalization constant of gluon mass. With the above definitions 
of the renormalization constants, on inserting Eq. (4.24) into Eqs. (4.11) and (4.22) , the ghost 
particle propagator and gluon propagator can be renormalized, respectively, in such a manner 

iA ab (k) = Z 3 iA a ^(k) (4.26) 

and + 
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where 



and 



iDll(k) = Z,iD R fJk) (4.27) 



.(/» , 1, ..afcr — k^k v /k 2 Z' 3 a R k^,k v jk 2 



are the renormalizcd propagators in which m R , ~p R and /i R are the renormalized masses, a R 
is the renormalized gauge parameter, £l R (k 2 ), IL]J(A: 2 ) and II^(fc 2 ) denote the finite corrections 
coming from the loop diagrams. They are defined as 

m R = z m lm > a R = Z 3 la , ~Pr = Mfl = \fz~3^, (4 30) 

n R {k 2 ) = z 3 fc(k 2 ), u T R {k 2 ) = z 3 [k 2 nt(k 2 ) + m 2 uuk 2 )}, n L R {k 2 ) = z' 3 w L {k 2 ). 

The finite corrections above are zero at the renormalization point v. As we see from Eq. 
(4.29), the longitudinal part of the gluon propagator, except for in the Landau gauge, needs 
to be renormalized and has an extra renormalization constant Z' 3 . This fact coincides with the 
general property of the massive vector boson propagator (see Ref. (8), Chap.V). From Eqs. 
(4.23)-(4.25) , it is easy to find that the longitudinal part in Eq. ( 4.22) can be renormalized as 

Z 3 a R [l + n R (k 2 )]A R (k 2 ) (4.31) 



fc 2 [i + n L (fc 2 )] -n 2 + is 

where 

a » ( ^ " f[i + n «wi-^ + fe <" 2) 

which appears in Eq. (4.28) and the renormalization constant Z' 3 can be expressed as 

^ = [l + 4ii T M]- 1 (4-33) 

If choosing v = /j, R , we have 

Z 3 = Z 3 (4.34) 



5. GAUGE-INDEPENDENCE AND UNITARITY 



This section serves to prove that the S-matrix elements evaluated by the massive gauge field 
theory are independent of the gauge parameter. That is to say, the gauge-dependent spurious 
pole appearing in the ghost particle propagator and the longitudinal part of the gauge boson 
propagator as shown in Eqs. (4.10) and (4.14) would not contribute to the S-matrix elements. 
This fact just ensures the unitarity of the S-matrix. According to the reduction formula' 11 !, 
the S-matrix elements can be obtained from the corresponding Green's functions. So, we first 
examine how the Green's functions are dependent on the gauge parameter. 

Let us start from the generating functional of Green's functions given in Eqs. (2.4) and (2.5). 
Since the fcrmion field in the generating functional is not related to the gauge parameter, for 
simplicity of statement, we will omit the fermion field functions in the generating functional 
and rewrite the functional in the form 

Z[J, K,K] = ± JV[A.C.C]exp{iS + i J d 4 x{-±.(d»A«) 2 + J-MJJ , , 

+K a C a + C a K a^ + ■ j d 4 xd 4 yC a (x)M ab (x , y)C b '(y)} 1 ' ; 
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where 

S = J d 4 x{-^F a ^F^ + ^m 2 A af *A a ^} (5.2) 

and 

M ab (x,y)=d»[Df(x)5 4 (x-y)} (5.3) 

in which V ab (x) was defined in Eqs. (2.6) and (2.7). 

When we make the following translation transformations in Eq. (5.1) 



C_ a {x) -» C a {x) - / d 4 y{M' 1 ) ab (x 1 y)K b (y) 
C a {x) -» C a (x) - / d 4 yK b (y)(M^ 1 ) ba (y, x) 



k„.ifbl„.\C\/r-\\ba( n . „\ (5-4) 



and complete the integration over the ghost field variables, Eq. (5.1) will be expressed as 

Z[J, K, K] = e -iI^yK^) { M-^ {x ,y,6/i6J)K\y) z[J] ^ 

where Z[J] is the generating functional without the external sources of ghost fields' 8 ' 12 !. 

Z[J] = 1 y 2>(,4)A F [4ezp{iS + i | d*a;[--L(^A-) 2 + J°»A<$} (5.6) 

in which 

A F [A] = detM[A] (5.7) 

here the matrix M[A] was defined in Eq. (5.3). From Eq. (5.5), we may obtain the ghost 
particle propagator in the presence of the external source J 



iA ab \x v J]- s zvxm 

= i(M-% b [x,y, 1 $ J ]Z[J] 



5K*(x)5K*(yyK=K=0 (5.£ 



The above result allows us to rewrite the W-T identity in Eq. (4.1) in terms of the generating 
functional Z[J] when completing the derivative with respect to u bu (y) and setting K a (x) = 

T?(y) = 0, 
where 

D b ?(y) = V b *(y)-^d,5 bd (5.10) 

is the ordinary covariant derivative. On completing the differentiations with respect to the 
source J, Eq. (5.9) reads 

j,fV[A}A F [A}exp{iS + ifd 4 x[-^Ap 2 + J^}} 

x[J d i y J b ^{y)D b ll c {y){M- 1 Y a {y,x) - l^A^x)] (5.11) 
= 

By making use of Eqs. (5.3), (5.8) and (5.10), the ghost equation shown in Eq. (4.2) may be 
written as 

M°°[x, lyKM-^tcy, ^j]Z[J] = 5 ab 8 4 (x y)Z[J] (5.12) 

When the source J is turned off, we get 

M ac (x)A cb (x -y) = S ab S 4 (x - y) (5.13) 
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This equation only affirms the fact that the ghost particle propagator is the inverse of the 
matrix M. 

Now we are in a position to describe the proof of the unitarity mentioned in the beginning 
of this section. To do this, it is suitable to use the generating functional written in Eq. (5.6) 
and the W-T identity shown in Eq. (5.11) because the S-matrix only has gluon external lines, 
without any ghost particle external line to appear. For simplifying statement of the proof, in 
the following, we use the matrix notation to represent the integrals. In this notation, Eqs. (5.6) 
and (5.11) are respectively represented as' 8 ' 

Z[J] F = ^J V(A)A F [A]e i ^-i F2+J - A ^ (5.14) 



and 



1 f V(A)A F [A]e^-i F2 + J - A Hj b D bc (M^) ca - -Ljr o ] = (5.15) 
N J y/a 

where we have defined 

F a ee ^A«(x) (5.16) 

with F corresponding to the gauge a, the subscript a, b or c stands for the color and/or the 
Lorentz indices and the space-time variable, and the repeated indices imply summation and/or 
integration. 

Let us consider the generating functional in the gauge a + Aa where Aa is taken to be 
infinitesimal 

Z[J] f+ af = ^J V(A)A F+AF [A}e^-^ F+AF ^ +J ^ (5.17) 



In the above, 



e 



-±(F+AF) 2 _ --fF' 



[1 + ^ 2 ] (5.18) 



A F+AF [A] =detM F+AF (5.19) 
According to the definitions given in Eqs. (5.3), (2.6) and (2.7), it is seen that 

M F b +AF = Mf + 6 ab Aam 2 (5.20) 

Therefore, 

A F+AF [A] = det[M F (l + Aam 2 M F 1 )} 

= detM F e Trln( > 1 + Aam2M F 1 ) ' (5.21) 
= A F [A\[\ + AarrfTrMp 1 ] 

Upon substituting Eqs. (5.18) and (5.21) into Eq. (5.17), we obtain 

Z f+ af[J] = iJV(A)A F [A]e^-i F2 +^ , , 

x{l+ l -^F 2 + Aam 2 TrM F 1 } ( ' 

For further derivation, it is necessary to employ the W-T identity described in Eq. (5.15). 
Acting on Eq. (5.15) with the operator ^Aaa~^F a [^j] and noticing 

jp Ke UcAc — iF f— 1— — e u <= A c — -$—F \JL] e iJcA c 



e »-A{iJ b F a [A] + »m} 



" J (5.23) 



we have 
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j ! JV(A)A F [A}e i ^-i F2 + JA y^{iJ b D bc [A](Mp 1 ) ca F a [A} 
+ %^ hc [A](M F 1 ) ca -^F 2 } = 

Adding Eq. (5.24) to Eq. (5.22) and considering 

^D bc (M^) ca = ^(dDUiM^U = ^[M - ^] ac {M^) ca 
= ^Tr[l - n^Mp 1 ] 

where Eqs. (5.10) and (5.3) has been used, one may reach the following result 

Z f+ af[J] = ± JV(A)A F [A]e« s W-? F2 + J - A }{l + zAS + U a [^D ab {M^) bc F c ]} 
= £ JV(A)A F [A]e^ s ^+ AS -i F2 + J - A '} 



(5.24) 



(5.25) 



(5.26) 



where 



A' a =A a + ■^^=D ab (Mp 1 ) bc F c (5.27) 



AS = - % -^Tr\l + l i 2 Mp 1 ] (5.28) 



in which 



TrMp 1 = J d 4 a;A aa (0) = const. (5.29) 

Since the AS is a constant (even though it is infinite), it may be taken out from the integral 
sign and put in the normalization constant N. Thus, Eq. (5.26) will finally be represented as 

Z f+ af[J] = j^J V(A)A F [A}e i ^-i F2+J - A '^ (5.30) 

In comparison of Eq. (5.30) with Eq. (5.14), it is clear to see that the difference between 
the both generating functionals merely comes from the vector potentials in the external source 
terms, while, the remaining terms belonging to the dynamical part in the both generating 
functionals are completely the same. This indicates that any change of the gauge parameter 
can only lead to different field functions in the source terms of the generating functional. 
According to the equivalence theorem, the different field functions in the source terms does 
not influence on the S-matrix elements, it can only affect the renormalization of external lines 
for the Green's functions and wave functions' 8-12 ' . This point will be explained more specifically 
in the following. 

The n-point gluon Green's functions computed from the generating functionals Z F [J] and 
Z f +af[J] are represented in the position space as 

G F {x u x 2 ,-- ;x n ) = (0 I T{A(x 1 )A(a; 2 ) • ■ • A(x n )} | 0) (5.31) 

G f+ af(xi,x 2 , ■ ■ ■, x n ) = (0 | T{A' ( Xl )A' (x 2 ) ■ ■ ■ A' (x n )} | 0) (5.32) 

where A(x^) and A'(ajj) denote the field operators corresponding the gauges F and F + AF 
and Xi designates the space-time for the i-th particle. Here and afterward, we adopt the matrix 
notation to represent the field functions and Green's functions, therefore, the Lorentz and color 
indices are suppressed. In light of the renormalization of the field operators' 11 ' 

A a ^{x) = z\a°£(x) (5.33) 
A'^(x) = Z F+AF A'^(x) (5.34) 
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where R marks the renormalized quantities and Z F and Zp+AF are the renomalization constants 
given in the gauges F and F + AF respectively, the renormalization of the above Green's 
functions, in the momentum space, can be written as 

G f (fci,fc 2 ,...,fc„) = ZjG$(k 1 ,k 2 ,...,k n ) (5.35) 



G F+AF (hk 2 , ■-, K) = z F+AF G F+AF (k 1 k 2} k n ) (5.36) 

It is well-known that the Green's functions computed from the generating functionals Z F [J] 
and Z F+AF [J] have different external lines, but the same internal structure. Therefore, by the 
equivalence theorem and noticing Eqs. (5.35) and (5.36), we have 

n 

II lim fe 2_ >m 2 (kf - m 2 R )G F+AF (k 1 k 2 , ...,k n ) 

: =1 n ™ ( 5 - 37 ) 

= z f+af/ z f II l ^k^ m %(ki - m 2 R )G F (k h k 2 , k n ) 

i=l 

The propagators given in the gauges F and F + AF are respectively renormalized in such a 
manner ^ 

D F (ki) = Z F D R (h) = D T F {ki)T + R F {h) (5.38) 



D F+AF (ki) = Z F+AF D R (k t ) = D T F+AF {ki)T + R F+AF (h) (5.39) 

where T denotes the transverse projector, D F (ki) and D F+AF (ki) come from the transverse 
parts of the propagators D F (ki) and D F+AF {ki) and have a physical pole at kf — m 2 R 

D F(h) = t^-t (5.40) 

K i m R 



D T F+AF (h) = ^f^r (5.41) 

while, R F (ki) and R F + AF (k,i) represent the remaining parts of the propagators D F {ki) and 
D F + AF {ki) which are regular at the physical pole, therefore, 



lim (kf - m%)R F (ki) = lim (kf - m R )R F+AF (ki) = (5.42) 



According to the reduction formula, the S-matrix elements for multi-gluon scattering which are 
evaluated in the gauges F and F + AF may be respectively represented as 



S F (h, k n ) = Z F 2 Yl lim fc 2^ m 2 A (ki)(kf - m R )G F (ki, fc„) 

1=1 ' (5-43) 
= Z 2 n lim fe 2^ ro 2 A (h)D];(h)- 1 G F (k 1 , k n ) 



i=i 



S F+AF (ki,...,k n ) = Z F 2 AF J] ^k 2 ^ m 2 R A o( k i)(kf - m 2 R )G F+AF (k 1 , ...,k n ) 

1=1 * (5-44) 
= z f+ A f U^ m k 2 ^ m % A o(ki)D^ +AF (ki) 1 G F+AF (k 1 ,... 1 k n ) 

where Ao(ki) is the free wave function of the i-th gluon which represents the state of transverse 
polarization and therefore is free of the gauge parameter. On substituting Eq. (5.37) into Eq. 
(5.44) and noticing Eq. (5.43), we find 
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SF+AF{h, fen) = S F {kl, k n ) 



(5.45) 



which shows that the S-matrix elements arc independent of the gauge parameter. The gauge- 
independence of the S-matrix elements implies nothing but the unitarity of the S-matrix because 
the gauge-dependent spurious pole which appears in the longitudinal part of the gluon propa- 
gator and the ghost particle propagator and represent the unphysical excitation of the massive 
gauge field in the intermediate states must eventually be cancelled out in the S-matrix ele- 
ments. From the construction of the theory, the cancellation seems to be natural. In fact, in 
the massive Yang-Mills Lagrangian, all the unphysical degrees of freedom have been restricted 
by the constraint conditions imposed on the gauge field and the gauge group. When these 
constraint conditions are incorporated in the Lagrangian, the theoretical principle we based on 
would automatically guarantee the cancellation of the unphysical excitations. This conclusion 
drawn from the above general proof can be checked by practical perturbative calculations, as 
will be demonstrated in a subsequent paper. 



6. REMARKS 



In the last section, we would like to make some remarks on the BRST- external source terms 
introduced in Eq. (3.4). Ordinarily, to guarantee the BRST-invariance of the source terms, 
the composite field functions A<I>i are required to have the nilpotency property SA&i = 
under the BRST-transformations' 8-11 !. For the massless gauge field theory, as one knows, the 
composite field functions are indeed nilpotent. This nilpotency property is still preserved for 
the massive gauge field theory established in the Landau gauge because in this gauge the BRST- 
transformations are the same as for the massless theory. However, for the massive gauge field 
theory set up in the general gauges, we find SA$i ^ 0, the nilpotency loses, since in these gauges 
the ghost field acquires a spurious mass fj,. In this case, as pointed out in section 2, to ensure 
the BRST-invariance of the source terms, we may simply require the sources u, to satisfy the 
condition denoted in Eq. (3.7). The definition in Eq. (3.7) for the sources is reasonable. Why 
say so? Firstly, we note that the original W-T identity formulated in Eq. (3.2) does not involve 
the BRST- sources. This identity is suitable to use in practical applications. Introduction of 
the BRST source terms in the generating functional is only for the purpose of representing the 
identity in Eq. (3.2) in a convenient form, namely, to represent the composite field functions 
in the identity in terms of the differentials of the generating functional with respect to the 
corresponding sources. For this purpose, we may start from the generating functional defined 
in Eq. (3.4) to re-derive the identity in Eq. (3.2). In doing this, it is necessary to require 
the source terms UiA$i to be BRST- invariant so as to make the derived identity coincide 
with that given in Eq. (3.2). How to ensure the source terms to be BRST-invariant? If the 
composite field functions A$j are nilpotent under the BRST-transformation, SA$i = 0, the 
BRST-invariance of the source terms is certainly guaranteed. Nevertheless, the nilpotency of 
the functions A<j>i is not a uniquely necessary condition to ensure the BRST- invariance of the 
source terms, particularly, in the case where the functions A$i are not nilpotent. In the latter 
case, considering that under the BRST- transformations, the functions A$j can be, in general, 
expressed as <5A$i = £<&i where the <J>i are some nonvanishing functions, we may alternatively 
require the sources Ui to satisfy the condition shown in Eq. (3.7) so as to guarantee the source 
terms to be BRST- invariant. Actually, this is a general trick to make the source terms to be 
BRST-invariant in spite of whether the functions A<I>i are nilpotent or not. As mentioned before, 
the sources themselves have no physical meaning. They 0X6, ClS ci mathematical tool, introduced 
into the generating functional just for performing the differentiations. For this purpose, only a 
certain algebraic and analytical properties of the sources are necessarily required. Particularly, 
In the differentiations, only the infinitesimal property of the sources are concerned. Therefore, 
the sources defined in Eq. (3.7) are mathematically suitable for the purpose of introducing 
them. The reasonability of the arguments stated above for the source terms is substantiated by 
the correctness of the W-T identities derived in section 4. Even though the identities in Eqs. 
(4.1) and (4.2) are derived from the W-T identity in Eq. (3.8) which is represented in terms of 
the differentials withe respect to the BRST-sources, they give rise to a correct relation between 
the gluon propagator and the ghost particle one as shown in Eq. (4.8). The correctness of 
the relation in Eq. (4.8) may easily be verified by the free propagators written in Eqs. (4.10) 
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and (4 14). These propagators were derived in paper I by employing the perturbation method, 
without concerning the BRST-source terms and the nilpotency of the BRST- transformations. 
A powerful argument of proving the correctness of the way of introducing the BRST-sources is 
that after completing the differentiations in Eq. (3.8) and setting the BRST-sources to vanish, 
we immediately obtain the W-T identity in Eq. (3.2) which is irrelevant to the BRST-sources. 
Therefore, all identities or relations derived from the W-T identity in Eq. (3.8) are completely 
the same as those derived from the identity in Eq. (3.2). An important example of showing this 
point will be presented in Appendix where the identity in Eq. (5.11) which was derived from 
the W-T identities in Eqs. (3.8) and used to prove the unitarity of the theory can equally be 
derived from the generating functional in Eq. (5.6) which does not involve the BRST-sources. 
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8. APPENDIX 

To confirm the correctness of the identity given in Eq. (5.11), we derive the identity newly 
by starting from the generating functional written in Eq. (5.6). The generating functional in 
Eq. (5.6) was directly derived from the massive Yang-Mills Lagrangian by the Faddeev-Popov 
method of quantization' 12 !. Let us make the ordinary gauge transformation 

SA a = D ab 0b (A1) 

to the generating functional in Eq. (5.6). Considering the gauge-invariance of the functional 
integral, the integration measure and the functional Aj?[A] = detM[A], we get' 8 ' 12 ' 

5Z[J] = ^jD(A)A F [A}Jd i y[J^(y)+m 2 A^(y) 
-^Ald^{y)e c {y) exp{iS + i J d^-^A^f + J a »A«]} (A2) 

= 

According to the well-known procedure, the group parameter 9 a {x) in Eq. (A2) may be deter- 
mined by the following equation' 5 ' 9 ! 

M ab (x)9 b (x) = d^Vf{x)e\x)) = A» (A3) 

where A a (x) is an arbitrary function. When setting A a (ir) = 0, Eq. (A3) will be reduced to the 
constraint condition on the gauge group (the ghost equation) which is used to determine the 
9 a (x) as a functional of the vector potential A a (x). However, when the constraint condition 
is incorporated into the action by the Lagrange undetermined multiplier method to give the 
ghost term in the generating functional, the 9 a (x) should be treated as arbitrary according to 
the spirit of Lagrange multiplier method. That is why we may use Eq. (A3) to determine the 
functions 9 a (x) in terms of the function X a (x) . From Eq. (A3), we solve 

9 a (x) = j d 4 x(M- 1 ) ab (x-y)X b (y) (A4) 

Upon substituting the above expression into Eq. (A2) and then taking derivative of Eq. (A2) 
with respect to A a (x), we obtain 

i / D(A)A F [A] J d 4 y[J b »(y) + m 2 A^(y) 
~a®yA b l/ (y)d£]D b l c (y)(M~ 1 ) ca (y - x) exp{*S+ 

i/dM-£(0Mj;) 2 + J»M«]} lAOj 
= 

According to the expression denoted in Eq. (2.7) and the identity f<> cd A°v A^ — o, it is easy to 
see 
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A b ^y)D b ;(y)(M- 1 r(y - x) = A b ^(y)d^(M~ 1 ) ba (y - x) (A6) 

By making use of the relation in Eq. (5.10), the definition in Eq. (5.3) and the equation in Eq. 
(5.12), we deduce 

^9yA b (y)d^D b ^(y)(M~ 1 ) ca (y - x) 
= ±d»A b v (y)6 4 (x -y)- m 2 d- y Al(y)(M- l ) ba {y - x) {A '> 

On inserting Eqs. (A6) and (A7) into Eq. (A5), we obtain an identity which is exactly identical 
to that given in Eq. (5.11) although in the above derivation, we started from the generating 
functional without containing the ghost field functions and the BRST-sources and , therefore, 
the derivation does not concern the nilpotency of the composite field functions appearing in the 
BRST-sourcc terms. This fact indicates that the W-T identities derived in section 3 are correct 
and hence the procedure of introducing the BRST-invariant source terms into the generating 
functional is completely reasonable. 
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